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摘 要 : 在 牧区 和 灌溉 农业 区 ,大 量 含 氮 的 畜 禽 排泄 物 和 氮肥 从 土壤 进入 地 表 水 和 地 下 水 ,是 流域 面 源 污染 的 主 
要 来 源 。 河 岸 潜流 带 是 削减 氛 素 污染 负荷 的 有 效 屏障 ,厘清 河岸 潜流 带 对 氮 素 的 迁移 转化 和 去 除 作 用 是 控制 流域 
氮 素 污染 的 关键 。 本 研究 选取 位 于 典型 草原 牧区 的 锡 林 河 上 游 河 段 ,开展 了 夏 汛期 河水 与 河岸 地 下 水 的 水 位 、 氨 
ib( NH; ) 和 硝酸 盐 ( NO; ) 浓 度 ,以 及 相关 环境 因子 的 连续 监测 ,并 利用 FEFLOW 建立 了 河岸 潜流 带 水 流 和 毛 素 溶 
质 反应 运 移 模 型 。 利 用 实测 数据 拟 合 的 模型 能 够 准确 再 现 河岸 潜流 带 水 位 和 两 种 主要 氮 素 浓度 的 动态 变化 。 结 
果 表 明 :(1) 夏 汛期 河岸 带 氮 素 污染 风险 较 高 ,河岸 带 NH 浓度 从 降雨 前 的 0.2 mg'I 升 高 到 降雨 后 的 7.23 mg L7, 
NO; 浓度 从 1 mg LFH A 8.27 mg*L"'。(2) 实测 和 模拟 结果 均 显 示 潜 流 带 中 氮 素 动态 与 降雨 .地表 水 -地 下 水 交换 
等 水 文 过 程 密切 相关 ,上 且 NH; I NO; 浓度 对 暴雨 事件 的 响应 机 制 不 同 ,(3) 降雨 期 间 ,流动 性 较 强 的 NO; 在 淋 滤 作 
下 从 河水 和 地 表 入 渗 进 入 河岸 带 , 导 致 浓度 显著 升 高 。 同 时 ,降雨 事件 加 强 了 河水 -地 下 水 的 交换 作用 ,通过 控 
制 营养 物质 的 输入 影响 氮 素 生物 地 球 化 学 循环 ,从 而 调节 河岸 潜流 带 NH; 和 NO; 浓度 的 变化 。 本 研究 初步 揭示 了 
牧区 河岸 带 对 于 氮 素 的 水 文 和 生物 地 球 化 学 过 程 的 缓冲 作用 机 制 , 为 牧区 氮 素 污染 控制 提供 了 科学 参考 。 

关键 词 : 河岸 潜流 带 ; ARKE; 水 文 过 程 ; 地 下 水 模拟 ; 锡 林 河 


氮 素 引起 的 富 营 养 化 是 全 球 范围 关注 的 水 环 。 污染 物 浓度 与 降雨 .流量 的 动态 关系 ,利用 滞后 现 
境 问 题 n"。 农 业 生 态 系 统 中 合成 肥料 的 过 量 施 用 和 ” 和 象 (hysteresis ) 解 释 氮 素 的 来 源 和 途径 , 即 氮 素 是 来 
畜牧 养殖 业 中 畜 禽 数 尿 的 大 量 排出 是 氮 素 面 源 污  ” 自 于 涨 水 期 通过 地 表 径 流 快速 进入 河流 或 是 在 落 
染 的 重要 来 源 。 大 气 ,地表 和 土壤 中 的 过 量 氮 素 在 。 ”水 期 由 地 下 水 排泄 至 河流 ,但 是 仅 依 靠 河水 的 流 
水 文 过 程 的 驱动 下 进入 水 体 , 从 而 造成 水 环境 污 ” 量 和 氮 素 浓度 监测 及 滞后 分 析 不 足以 准确 识别 氮 
W ,硝酸 盐 ( NO; ) 和 氨 盐 ( NH; ) 是 水 体 中 最 主要 的 。 素 的 迁移 转化 ,需要 结合 氮 素 及 相关 环境 因子 在 河 
氮 素 形态 "。 河 岸 潜流 带 作 为 地 表 水 和 地 下 水 的 混 。” 岸 带 的 分 布 进一步 明晰 其 周转 过 程 。 
合 区 域 ,其 中 的 水 量 交换 促进 了 物质 和 能 量 的 交 大 多 数 研 究 表明 ,河岸 潜流 带 是 硝酸 盐 的 汇 ， 
换 , 并 形成 了 明显 的 氧化 还 原 梯度 和 多 样 的 生物 群 ”其 中 发 生 的 硝化 ` 反 硝化 . 厌 氧 氮 氧 化 等 过 程 取决 
沙 , 因 此 被 普遍 认为 是 营养 物质 转化 ,特别 是 氮 素 ”于 潜流 融 中 的 氮 素 供应 水 化 学 特征 以 及 水 文 条 件 
循环 的 热点 区 域 59。 厘 清 两 种 氮 素 在 河岸 潜流 带 “(例如 洪水 特征 、 地 下 水 位 ) 等 '"。 目 前 ,研究 者 对 
的 迁移 转化 和 来 源 是 氮 素 污染 负荷 控制 的 关键 。 灌溉 农业 区 附近 的 河流 潜流 带 开 展 了 广泛 研究 , 例 

有 研究 表明 ,夏季 汛期 是 氮 素 转化 的 热点 时 如 ,Darwiche-Criado 47 RIM, NO; 在 丰 水 期 与 枯 
期 ,降雨 径流 事件 是 流域 氮 素 面 源 污染 的 主要 驱动 水 期 的 来 源 不 同 , 且 对 水 文 过 程 的 响应 不 同 。 潘 俊 
因子 ,刺激 了 和 氮 循 环 的 周转 ”。 相 关 研 究 通 过 分 析 。 ”等 "5" 发现 了 潜流 带 的 氧化 还 原 条 件 影响 了 氮 素 的 
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生物 地 球 转化 过 程 , 反 硝化 、 大 氧 氨 氧化 转化 强度 
在 一 定 范围 内 存在 空间 差异 。 但 是 ,目前 国内 人 研究 
对 牧区 河岸 之 的 氮 素 迁移 转化 规律 认识 仍 有 不 足 ， 
大 多 人 研究 仅 关注 NO; 在 反 硝 化 作用 下 的 迁移 转化 
模拟 “”。 基 于 此 ,本 研究 基于 汛期 高 频率 水 位 ,水 
体 氮 素 浓度 及 其 相关 环境 因子 的 监测 ,采用 模拟 软 
件 FEFLOW2 构建 河岸 潜流 带 水 流 - 氮 素 迁移 转化 
模型 。 忆 在 分 析 河 流 及 河岸 潜流 带 中 氮 素 动态 及 
其 对 水 文 过 程 的 响应 ,探究 河岸 潜流 融 中 两 种 主要 
AUR TRES NO, , NH; 的 来 源 和 运 移 机 制 , 以 期 为 河 
岸 带 管 理 及 流域 氮 素 污染 治理 提供 借鉴 与 帮助 。 


1 材料 与 方法 


1.1 研究 区 概况 

人 研究 区 域 位 于 内 蒙古 高 原 东部 的 锡林郭勒 草 
原 , 地 处 我 国 草原 从 东部 半 湿 润 草 多 草 原 区 向 西北 
半 干 旱 典型 草原 区 过 渡 地 带 ( 图 1)。 锡 林 河 是 锡 林 
郭 勒 草原 的 母亲 河 , 全 长 268 km, 从 东南 向 西北 流 
动 。 锡 林 河 流域 面积 10542 km’, 年 降水 量 为 250~ 
400 mm, 全 年 最 大 降水 量 出 现在 夏季 汛期 的 7 一 8 
月 ,多 年 平均 蒸发 量 1900 mm", He PAKS a 
孔隙 水 为 主 , 富 水 性 较 强 ,上 层 潜水 含水 层 以 上 更 
新 统 的 细 砂 为 主 。 锡 林 郭 勒 草原 是 我 国 北方 重 
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要 的 牧场 。 
12 原 位 定点 监测 试验 

于 2020 年 夏季 汛期 ,在 位 于 中 国 科学 院内 蒙古 
草原 生态 系统 定位 研究 站 附近 的 锡 林 河 上 游 河 段 
开展 野外 定点 监测 试验 。 监 测 时 段 为 6 月 15 日 、7 
月 2 日 至 8 月 14 日 .9 月 20 日 。 人 研究 河 段 总 长 110 
m, 河 宽 平 均 8.5 m。 监 测 期 内 降雨 量 共 180.3 mm, 
其 中 两 次 较 大 降雨 ,分 别 为 7 月 11 一 13 日 连续 3 d 
降雨 共 47 mm 和 8 月 8 日 的 降雨 64 mm。 监 测 期 内 
河水 深 0.3~0.7 m, 河 岸 地 下 水 埋 深 为 0.5~1.1 m。 分 
别 在 河道 及 河岸 共 布 设 9 个 监测 井 ( 图 1), 河 道内 间 
隔 5 m 均 匀 布 设 3 个 监测 井 (C1~C3), 河 道 左岸 (1L~ 
3L) 和 右岸 (1R~3R ) 分 别 布设 3 个 监测 井 , 位 于 距离 
河道 7~10 m 位 置 。 河 岸 监测 管 深度 为 地 表 以 下 2 
m, 河 道 监测 管 深 度 为 河床 表面 以 下 1 m。 选 用 直径 
为 50 mm 的 PVC 管 作为 地 下 水 监测 井 管 ,在 距离 管 
底 10~15 cm 的 位 置 均匀 钻 孔 ,并 用 400 目 尼 龙 网 包 
庄 , 防 止 洲 泥 堵塞 。 

监测 前 期 ,采用 全 站 仪 测 量 研 究 区 地 形 。 监 测 
期 间 , 采 用 水 位 计 及 水 尺 逐 日 监测 地 下 水 位 及 河水 
位 。 使 用 流速 仪 测量 河水 流速 ,根据 河 宽 换 算 流 
A. f S~7 d 采 用 便携 式 水 质 分 析 仪 (Multi 3630 
IDS SETG ) 监 测 河水 和 地 下 水 的 温度 溶解 氧 (Dis- 
solved Oxygen, DO) .pH 等 水 质 参数 。 同 时 采集 河 
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图 1 研究 区 及 研究 河 段 监 测 井 布设 


Fig. 1 Overview of the study area and layout of the monitoring wells in the studied river section 
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水 和 地 下 水 水 样 , 置 入 100 mL 棕色 采样 瓶 AH. 
避 光 保存 ,防止 发 生生 物化 学 转化 。 将 采集 水 样 在 
实验 室 通过 直径 0.4 um 滤 膜 过 滤 ,根据 生态 环境 部 
推荐 的 测定 方法 ,使 用 Vario TOC 分 析 仪 测定 溶解 有 
机 碳 (Dissolved Organic Carbon, DOC)TKHE^ , (EH 
紫外 分 光 光 度 计 (SHIMADZU UV-2600) 测 定 无 机 
AU BE : NH; 采用 纳 氏 试剂 分 光 光 度 法 测定 2 ; 
NO; 采用 紫外 分 光 光 度 法 测定 1。 
1.3 河岸 潜流 带 水 流 和 溶质 运 移 模 型 构建 

人 研究 区 河岸 潜流 带 通 过 降水 和 河流 渗透 得 到 
补给 ,河水 和 地 下 水 的 交换 取决 于 河流 与 地 下 水 位 
之 间 的 水 力 梯度 和 沉积 物 的 渗透 性 能 。 放 牧牛 群 
排泄 的 业 尿 是 研究 区 主要 的 氮 素 来 源 , 在 水 文 和 生 
物 地 球 化 学 的 共同 作用 下 可 能 发 生 硝化 ` 反 硝化 等 
反应 过 程 。 本 研究 基于 FEFLOW 建立 了 人 研究 河 段 
河岸 潜流 带 的 三 维 水 流 和 毛 素 溶质 反应 运 移 数 值 
模型 (图 2) ,根据 潜流 带 内 水 流 过程 和 毛 素 迁移 转 
化 过 程 ,采用 2020 年 夏季 汛期 的 地 下 水 位 和 地 下 水 
NH; , NO; 浓度 进行 模型 拟 合 。 模 型 长 65 m, 宽 40 
m, 深 度 4m, 采 用 三 角形 有 限 单元 网 格 进行 离散 ,在 
河道 区 域 加 密 。 根 据 实测 地 表 高 程 设置 模型 ,垂直 
分 为 4 层 , 每 层 由 1069 个 单元 格 构成 。 系 统 自 动 生 
成 每 天 9 个 模拟 时 间 步 长 ,共计 44 do 


高 程 /m 

m 1180.6 
m 1180.2 
= 1179.7 
@ 1179.3 
m 1178.9 
m 1178.5 
W 1178.1 
Wm 1177.6 
m 1177.2 
m 1176.8 


图 2 FEFLOW 中 的 数字 高 程 模型 
Fig. 2 Digital elevation model in FEFLOW 


河岸 潜流 带 中 水 流 采用 非 稳 定 流 潜水 方程 
MR: 
ð oH ð oH 
+ 
Ox (Kh Ox ) oy (KA oy ) 


à aH aH 
+ +o =u% 
人 


(1) 
式 中 :KK 天 是 沿 坐标 轴 x yz 方向 的 渗透 系数 
(m: d^) ,假设 坐标 轴 与 渗透 系数 的 方向 一 致 ;为 
地 下 水 位 (my) ;为 潜水 层 厚度 (m) it Jg SERI Cd) 5 o 
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6 期 酝 栋 元 等 :牧区 河岸 潜流 带 硝 酸 盐 气 和 毛毛 浓度 对 水 文 过 程 的 响应 机 制 939 


为 单位 时 间 水 平面 积 上 的 垂 向 补给 量 (m d) ;为 
给 水 度 。 
溶质 反应 运 移 模 型 采用 以 下 对 流 弥 散 方 程 
摘 述 2 : 
oe =V(DV,)-V(vc) +R, (2) 


SUF ie HK AYA J AYR BE Ing L') St AT Ta Cd); 
万 为 水 动力 弥散 系数 (ms 站 ,分 为 纵向 和 横向 , 
广 表 示 纵 向 弥散 系数 (m) ,六 表示 横向 弥散 度 
(m) ,可 假设 其 与 溶质 浓度 无 关 ;v 为 平均 地 下 水 流 
速 (m*s 1!);V 为 梯度 ;V, 为 浓度 梯度 ;RR 为 溶质 源 汇 
Ji(mg*m^:s'). 
D=av +D, (3) 
式 中 ;a 为 弥散 度 (m) ,分 为 纵向 和 横向 ,a 表示 纵 问 
弥散 度 (m) ,ar 表 示 横 问 弥 散 度 (m) ;Do 为 分 子 扩散 
Aw). 
NO; 移动 性 较 强 , 易 受 到 淋 深 作用 和 反 硝 化 影 
响 ,而 NH; 移动 性 较 低 , 易 发 生 吸 附和 硝化 等 作 
JH. AES) NH; 的 主要 来 源 是 有 机 毛 的 矿 化 ,而 研 
究 区 有 机 氮 来 源 丰 富 ,并 且 模 型 为 河岸 带 尺 度 ,不 
考虑 因 氮 氮 被 吸附 而 引起 的 迁移 范围 受 限 。 考 虑 到 
锡 林 河 上 游 河 兰 带 土壤 腐殖质 丰富 , 且 相 关 人 研究 表 
明 河岸 带 中 参与 氮 素 转化 的 微生物 种 类 丰富 *”， 
特别 在 温度 适宜 的 夏季 较为 活跃 。 因 此 ,本 研究 认 
为 微生物 的 生长 -死亡 过 程 并 不 是 影响 潜流 带 中 氛 
素 迁 移 转化 的 主要 因素 。 在 模型 中 将 微生物 作为 
有 充足 .稳定 来 源 的 量 ,并 未 考虑 微生物 的 生长 - 死 
亡 过 程 。 采 用 一 级 反应 动力 学 描述 公式 (2) 中 NH; 
BONO, 的 衰减 ”, 源 汇 项 如 下 : 
Ra =N KX Co (4) 
Ra =N nia EXC. (5) 
式 中 :Ri 为 NO; 的 源 汇 项 ;和 Vi 表示 土壤 淋 滤 进 入 地 
下 水 的 NO; 量 (mg.m”…s');K 表 示 由 各 类 反应 ,如 
有 反 硝化 导致 NO; 浓度 降低 的 衰减 系数 (d ); C、 为 
NO; 浓度 ;Rj 为 NH; 的 源 汇 项 ;Ni 表示 有 机 毛 矿 化 
Œ (mg m™: s~!) ;有 表示 由 吸附 、 硝 化 等 过 程 导致 
NH; 浓度 降低 的 衰减 系数 (d 7) 5 C... 为 NH 浓度 。 
在 研究 区 域内 的 河流 位 置 设 为 第 一 类 边界 
(Dirichlet 边 界 ) ,根据 逐日 监测 的 河水 位 确定 。 降 
水 量 作 为 源 汇 项 ,降水 入 渗 系 数 根据 区 域 包 气 带 下 
垫 面条 件 作为 经 验 值 计 算 。 课 题 组 前 期 研究 结 
表明 3, 人 研究 河 段 侧 向 潜流 交换 边界 位 于 距离 河道 
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10-15 mm 的 位 置 ,该 模型 前 后 边界 距 河流 约 为 20 m, 
此 处 河岸 地 下 水 与 河流 之 间 的 交换 非常 微弱 ,因此 
设置 为 零 通 量 边界 。 左 右边 界 可 能 有 地 下 水 流入 
流出 ,设置 为 流量 边界 。 垂 问 潜 流 交换 边界 位 于 河 
床 表面 以 下 1m 左右 ,该 模型 下 边界 位 于 河床 
表面 以 下 4m, 基 本 无 垂 向 交换 ,因此 模型 底部 设 为 
零 通 量 边界 。 河 水 NO; 及 NH; 浓度 分 别 作为 溶质 
模块 河流 的 定 浓度 边界 条 件 , 其 余 边界 的 浓度 均 为 
0。 以 2020 年 7 月 3 日 的 地 下 水 位 和 溶质 浓度 的 空 
间 分 布 作为 初始 条 件 ,根据 监测 期 内 实测 水 位 及 溶 
质 浓度 对 模型 进行 率 定 。 

运用 抽水 法 测定 河床 和 河岸 带 沉积 物 的 渗透 
系数 ” ,模型 参数 汪 ” 见 表 1。 

为 了 验证 模型 模拟 精度 ,使 用 均 方 根 误差 
(Root Mean Square Error, RMSE ) 和 拟 合 优 度 (Coeffi- 
cient of Determination , ) 作 为 评价 指标 ,公式 如 下 : 


(6) 


gai 2407-0 


5 Y y» 


(7) 


式 中 ;y 表 示 实 际 测量 值 ; y 表示 实际 测量 平均 值 ;y" 
表示 模型 模拟 值 ;i=1,2,… ,n,n 表示 数据 个 数 。 
1.4 模型 参数 敏感 性 分 析 

根据 模型 模拟 污染 物 溶质 特征 ,选取 了 5 种 对 
河岸 潜流 带 溶质 运 移 存 在 影响 的 参数 进行 敏感 性 
分 析 。 参 数 包括 : 纵 问 弥 散 度 a、 横 向 弥散 度 ar、 分 
子 扩散 系数 DD 一 级 反应 常数 KK,。 局 部 灵敏 度 分 
析 是 用 来 分 析 单 个 参数 变化 对 数值 模型 输出 结 
的 影响 程度 ,不 考虑 参数 之 间 的 耦合 作用 。 对 一 种 
参数 进行 灵敏 度 分 析 时 ,保持 其 他 参数 不 变 ,将 其 
取 值 分 别 增 大 10% 20% F40% ,分 别 记录 不 同 取 
值 情况 下 的 模拟 浓度 相对 于 调试 验证 后 模型 结 
(基本 情景 ) 的 变化 。 
1.5 河岸 潜流 带 NO; 去 除 效果 评价 

采用 RR 反映 河岸 潜流 带 内 NOS 的 产生 和 去 除 


as 


C 
Cs 
式 中 :CG 表示 河岸 潜流 带 中 NO0; 浓度 (mg: L’); C. 
表示 河水 中 NO; 浓度 (mg:L')。 
基于 Kunz 等 1 的 研究 方法 ,可 计算 河岸 潜流 带 
NO; BRR. HAASE : 
C-C 


Rcs (8) 


R=" n 100% (9) 
2 结果 与 分 析 
2.1 水 动力 学 过 程 


监测 期 内 发 生 两 次 较 大 规模 的 降雨 事件 :7 月 
11—13 日 连续 3 d MR 47 mm,8 月 8 日 降雨 64 mm 
(图 3)。 受 降雨 影响 ,河水 位 与 河岸 地 下 水 位 在 监 
测 期 内 频繁 波动 ,并 且 呈 现 相 似 的 变化 趋势 ,水 位 
变化 幅度 有 所 不 同 。7 月 9 日 河水 位 与 地 下 水 位 小 
幅 升 高 ,监测 区 内 无 较 大 降水 ,推测 可 能 是 由 于 上 
游 来 水 导致 。 在 7 月 11 一 13 日 连续 降水 期 间 , 地 下 
水 位 达到 监测 期 第 一 次 较 大 峰值 ,上 涨 近 0.30m , 河 
水 位 略 高 于 地 下 水 位 ,上 涨 0.45 m。 降 水 结束 后 ,水 
位 恢复 至 降雨 前 7 月 10 日 的 水 平 。8 月 8 日 第 二 次 
降雨 期 间 ,水 位 较为 平稳 , 仅 有 小 幅度 升 高 。 降 雨 
分 别 导 致 地 下 水 位 及 河水 位 上 涨 0.20 m,0.27 mo 
模拟 结果 表明 模型 能 够 较 准 确 得 描述 地 下 水 位 的 
实际 波动 过 程 (图 3) ,模拟 和 实测 的 地 下 水 位 拟 合 
程度 较 好 ,RMSE 为 0.04~0.07 m, R E 0.57~0.85 之 
间 。 模 拟 值 比 实测 值 相 比 略 有 出 和 人 ,可 能 是 由 于 不 
同 规模 降雨 时 的 河岸 入 渗 能 力 有 所 差异 ,导致 地 下 
水 位 的 涨幅 稍 有 不 同 。 

监测 期 内 频繁 的 水 位 波动 导致 了 河水 与 河岸 
地 下 水 水 力 梯 度 的 变化 。 无 降雨 时 水 力 梯 度 为 较 
小 的 负 值 ,表现 为 地 下 水 向 河流 微弱 的 排泄 趋势 
(图 4)。 当 降雨 发 生 时 ,河水 位 上 涨幅 度 较 大 ,导致 
河水 位 高 于 河岸 地 下 水 位 ,水 力 梯 度 发 生 道 转 , 表 
现 为 河水 入 渗 补 给 地 下 水 ,7 月 12 日 达到 监测 期 内 


表 1 模型 参数 取 值 


Tab.1 Model parameter values 


模型 参数 给 水 度 渗透 系数 纵向 
/(m? +b’) [m-a 弥散 度 /m 

参数 值 0.1 1.75~3.25 10 

取 值 依据 [35] 实测 [36] 


横向 分 子 扩散 一 级 反应 常数 /d" 
弥散 度 /m 系数 /ms ) NO; NH; 
1 1x10? 0.79 3.81 
[36] [37] [38-39] [38-39] 
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一 模拟 值 ”一 河水 位 


o 1180.2 
40 1179.9 
80 1179.6 
120 11793 


160 1179.0 


941 


降水 量 /mm 


1180.2 


6 期 醇 栋 元 等 :牧区 河岸 潜流 带 硝 酸 盐 
: 实测 值 
1180.2 1180.2 
3 1179.9 1179.9 
ii 1179.6 1179.6 LA 
1179.3 1179.3 
1179.0 1179.0 
1180.2 1180.2 
g 1179.9 1179.9 
x 
E 1179.6 


1179.3 
1179.0 + 


1179.3 
1179.0 


bud 1179.6 上 


40 1179.9 
80 1179.6 
120 11793 

1179.0 


图 3 降水 量 、 河 水 位 及 河岸 监测 井 实测 水 位 与 模拟 水 位 拟 合 


Fig. 3 Rainfall, river stage and fitting of measured and simulated water levels in river bank monitoring wells 
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图 4 监测 期 内 河水 与 河岸 地 下 水 之 间 水 力 梯度 随时 间 变 化 


Fig.4 Change of hydraulic gradient between river water and river bank groundwater during the monitoring period 


最 大 值 0.047。 降 雨 结 束 后 ,河水 位 下 降 较 快 ,恢复 


至 降雨 前 水 位 ,水力 梯 度 为 负 值 , 且 接近 0。 
2.2 氮 素 浓度 变化 


在 监测 期 内 ,河水 NO; 浓度 介 于 1.27~5.36 mg- 
L ,河岸 地 下 水 NO; 浓度 范围 为 0.29~8.27 mg:L'。 
河水 与 地 下 水 NO; 浓度 具有 相似 的 变化 趋势 ,呈现 


不 同 程度 的 波动 ,在 8 月 8 日 均 达到 峰值 ,除了 2R 点 
之 外 ,河水 浓度 略 低 于 地 下 水 浓度 (图 5)。7 月 12 
日 河岸 各 监测 井 NO; 平 均 浓度 由 降雨 前 的 0.93 mg， 
L EFT 228 mg*L', 随 后 下 降 至 1.47 mg-L';8 H 
8 日 经 历 64 mm 的 暴雨 后 ,河岸 监测 井 NO; 浓度 急 


剧 升 高 至 8.27 mg- L'o ARAM 


者 果 表明 模型 能 够 
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较 准 确 刻 画 河 岸 地 下 水 中 NO; 浓度 变化 趋势 ,模拟 
与 实测 结果 吻合 程度 较 高 ,RMSE 为 0.66~0.98 mg- 
L',R'fE0.64-0.89 之 间 。 模 型 略 低 估 了 1R、2R 和 
并 监测 井 在 8 月 8 日 的 NO; 浓度 峰值 ,而 在 3R 和 1L 
处 略 高 估 了 NO; 浓度 峰值 。 可 能 由 于 河岸 带 不 同 
位 置 的 氮 素 污染 负荷 输入 存在 差异 ,导致 实测 浓度 
较 模 拟 值 偏 高 或 偏 低 的 情况 出 现 。 


在 整个 监测 期 间 ,河水 NH; 浓度 介 于 0.15~ 
1.53 mg: LZ [E] ,河岸 地 下 水 NH} 浓度 范围 为 0.20~ 
7.23 mg*L'( 16). 97K NH; 浓度 处 于 较 低 水 平 ， 
趋势 平缓 , 主要 是 因为 河水 中 DO 充足 , NH; 易 发 生 
硝化 反应 转化 为 NO; *, 5; NO; 浓度 动态 相 比 , 河 
岸 监测 井 NH; 浓度 具有 不 同 的 变化 特征 。 不 同 监测 
井中 地 下 水 NH; 浓度 均 在 7 月 11 一 13 日 的 降雨 过 程 


。 实测 值 一 一 模 测 值 ”一 一 河水 浓度 
10 2R 10 [ 3R 
T 8 : 8 上 
El 6 6 + 
5 2 2+ . 2 
g . 
| E 0 L 
3 
E: 
E 
E 
D 
o 
Z 
图 5 不 同 监测 井 NO; 浓度 的 实测 值 与 模拟 值 拟 合 
Fig. 5 Fitting of measured and simulated NO; concentrations in different monitoring wells 
实测 值 — 模 测 值 —— 河水 浓度 
8 [ IR 8[ 2R 8[ 3R 
S s 
oD 
E 4 
EN 
É 2 
0 
8 
ty 6 
bp 
E 4 
EN 
X2 
Z o 


图 6 不 同 监测 井 NH; 浓度 实测 值 与 模拟 值 拟 合 


Fig.6 Fitting of measured and simulated NH; concentrations in different monitoring wells 
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中 明显 升 高 ,并 达到 峰值 (2.70~7.23 mg L^) ,降雨 结 
束 后 地 下 水 NH; 浓度 缓慢 下 降 到 监测 期 最 低 浓度 
(0.20~0.87 mg*L')。8 月 8 日 的 降雨 导致 地 下 水 NH; 
浓度 又 有 小 幅 上 升 趋 势 (1.48~2.29 mg-L), NH; TK 
度 模拟 值 变化 趋势 与 实测 结果 基本 吻合 ,RMSE 为 
0.61~0.70 mg: L'', R J 0.61-0.74,. 3L 监测 井 处 的 
NH; 模拟 浓度 小 于 实测 浓度 ,可 能 是 因为 放牧 导致 
的 局 部 高 浓度 氮 素 污染 物 排放 。 

由 敏感 性 分 析 结 果 ( 表 2) 可 知 , 纵 向 弥散 度 o 和 
横向 弥散 度 ar 参 数 取 值 增加 40% 时 ,引起 浓度 变化 
的 最 大 值 仅 为 0.067 mg L”, MATI RRA D HUE 
增加 40% 时 ,结果 变化 仍 小 于 0.001 mg- L'o EIEL 
上 3 个 参数 对 于 NO; 及 NH; 两 种 溶质 的 浓度 影响 较 
小 ,参数 灵敏 度 较 低 。 相 对 而 言 , NO; 的 衰减 系数 K, 
与 NH; 的 衰减 系数 已 参数 灵敏 度 较 高 。 当 玉 增 加 
40% 时 , NO; 峰值 浓度 降低 了 2.64 mg: L7; 25 K-I hn 
40% 时 , NH; 峰值 浓度 降低 了 1.13 mg:L'。 由 此 可 
见 , 反 应 参数 对 模型 结果 影响 较 大 , 即 河岸 带 的 生物 
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地 球 化 学 反应 对 于 氮 素 归宿 具有 重要 作用 。 
2.3 相关 水 环境 指标 变化 

在 监测 期 内 ,河水 pH 范围 为 7.78~8.31, 河 岸 地 
下 水 pH 范围 为 7.52~7.81, 均 为 弱 碱 性 。 河 水 DOC 
浓度 介 于 8.01~47.3 mg .二 ,河岸 监测 井 DOC 浓度 介 
于 1.83~29.72 mg L'o WE WJF DOC Yk REIME h 
现在 7 月 12 H , DOC Ye BE h KEk Bii Bj 11.21~25.18 
mg: L'HE 17.27~29.72 mg*:L'( 图 7)。 降 雨 结束 
后 ,DOC 浓 度 呈 现下 降 趋势 ,在 8 月 2 日 降 至 监测 期 
浓度 最 低 值 1.37~11.36 mgL-:。 河 水 中 DO 充足 
(7.74+0.34 mg L`) ,可 为 微生物 的 好 氧 呼吸 等 提供 
有 利 条 件 。 河 岸 地 下 水 中 的 DO 浓度 变化 范围 为 
0.2~1.7 mg: L, Ab T de AREE , CUBE JI Efe 3 
DOC 相似 ,同样 在 7 月 12 日 监测 到 最 大 值 。 


3 it it 


从 年 内 时 间 尺 度 来 看 ,夏季 汛期 河岸 潜流 带 的 


表 2 模型 参数 敏感 性 分 析 结 果 


Tab.2 Parameter sensitivity analysis results of the established model 


NO; 浓度 相对 于 基本 情景 的 减 小 量 /mg 


NH; 浓度 相对 于 基本 情景 的 减 小 量 /mg'L-) 


参数 增 量 
10% 20% 40% 10% 20% 40% 
纵向 弥散 度 om 0.011 0.029 0.067 0.007 0.013 0.035 
横向 弥散 度 or 0.001 0.003 0.007 0.001 0.001 0.002 
分 子 扩散 系数 D。 <0.001 <0.001 <0.001 <0.001 <0.001 <0.001 
NO; 衰减 系数 kK 1.029 1.410 2.642 - - - 
NH; $e AM Ko - - - 0.467 0.683 1.131 
一 一 DOC 浓 度 ” 一 一 DO 浓度 
FERES 2.5 307 2R 2.5 30 p 3p 2.5 
à 24 20 24 20 24 20 7, 
È 18 1.5 18 1.5 18 15 E 
R 12 10 12 10 12 10 x 
z 6 05 6 05 6 05 $ 
0 0.0 0 00 oh 1 1 1 L4 0.0 
ea PK FOI CP eat uu 
30 F iL 2.5 30 r 2L 2.5 30 2.5 
T 24 20 24 20 24 20 人 
È 18 1.5 18 1.5 18 1.5 E 
起 12 1.0 12 1.0 12 1.0 E: 
e 6 05 6 05 6 os 8 
0 0.0 0 0.0 0 0.0 
Fol e "i "i Fo Fd "i "i 
日 期 /月 -日 日 期 /月 -日 


图 7 监测 期 河岸 监测 井中 DOC 和 DO 浓度 变化 


Fig.7 DOC and DO concentration changes in riverbank monitoring wells during the monitoring period 
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氮 素 峰值 浓度 ( NO; 为 8.27 mg-L"', NH; 为 7.23 mg- 
Lg T UU Bk E ( NO; 为 5.07 mg- L’, NH; 为 
2.56 mg.I) 和 汛期 后 浓度 ( NO; 29 4.67 mg*L'、 
NH; 42.54 mg*L')。 主 要 是 因为 研究 区 全 年 降水 
量 少 , 且 集 中 于 夏季 汛期 ,地 表 和 土壤 中 积累 的 大 
量 NO; 在 降雨 的 冲刷 作用 下 进入 水 体 ” ,导致 河水 
和 洪流 带 中 NO; 浓度 显著 升 高 。 说 明 夏 汛 期 牧区 
河岸 带 氮 素 污染 的 风险 较 高 ,是 氮 素 污染 负荷 控制 
的 关键 时 期 % 5 ,但 是 NO 和 NH: 在 河岸 潜流 带 的 
迁移 转化 特征 不 尽 相 同 。 
3.1 汛期 河岸 潜流 带 中 NO; 的 迁移 转化 

试验 和 模拟 结果 表明 ,水文 过 程 对 河岸 潜流 带 
及 河流 中 氮 素 浓度 的 变化 起 了 主导 作用 ,但 是 控制 
NO; 和 NH; 迁移 转化 的 机 制 不 同 。 河 岸 潜流 带 
NO; 浓度 在 连续 3 d 强 度 较 小 的 降雨 事件 后 ,小幅 
上 升 而 后 下 降 , 可 能 是 土壤 表层 中 的 NO; 首先 被 淋 
滤 出 来 ”“ ,在 落水 期 浓度 缓慢 下 降 ”。 此 外 ,河岸 
潜流 带 DO 浓 度 升 高 为 硝化 反应 提供 了 条 件 ( 图 7)， 
促使 了 NH; 向 NO; 的 转化 “ ,导致 一 定 程 度 河岸 湾 
流 带 NO; 浓度 的 升 高 

在 8 月 8 日 的 暴雨 事件 后 ,河岸 潜流 带 NO; TK 
度 急 剧 升 高 至 监测 期 内 的 最 大 值 (图 5), 主 要 由 于 
强 降 雨 的 冲刷 作用 导致 土壤 上 层 中 累积 的 NOS 被 
淋 滤 出 来 ,一 部 分 随地 表 径 流 进入 河水 , 另 一 部 分 
可 能 直接 从 地 表 入 渗 进 入 河岸 地 下 水 。 由 于 NO; 
的 流动 性 较 强 ,降雨 冲刷 是 控制 河岸 带 NO; 浓度 的 
主要 作用 , 雨 强 的 大 小 决定 了 冲刷 作用 的 强 弱 以 及 
河水 补给 地 下 水 的 强 弱 *3。 
32 汛期 河岸 潜流 带 中 NH; 的 迁移 转化 

监测 期 内 的 河水 NH; 浓度 均 小 于 河岸 潜流 带 
地 下 水 浓度 。 在 7 月 11 一 13 日 连续 降雨 的 影响 下 ， 
河岸 潜流 带 NH; 浓度 与 DO .DOC 浓 度 一 同 逐 渐 升 
高 (图 6 和 图 7) ,于 7 月 13 日 达到 峰值 ,此 时 以 河水 
入 渗 补 给 河岸 地 下 水 为 主 , 河 水 与 地 下 水 的 水 力 梯 
度 也 达到 监测 期 内 最 大 (图 4)。 河 岸 潜流 带 中 NH; 
浓度 的 升 高 可 能 是 由 于 降雨 过 程 中 DO 和 DOC 随 河 
水 入 渗 进 入 地 下 水 ,促进 牛 群 排泄 物 中 的 有 机 毛 矿 
化 生成 NH; 。 研 究 表明 ,沉积 物 中 埋藏 的 有 机 碳 不 
能 很 好 地 促进 NH 的 释放 ,而 外 源 有 机 碳 作 为 重要 
的 电子 供 体 ,更 大 程度 地 促进 了 微生物 有 氧 呼吸 、 
有 机 毛 矿 化 等 过 程 ,这 可 能 是 降雨 期 间 潜流 带 中 


NH; 浓度 升 高 的 关键 驱动 因素 ”"。 虽 然 异 化 硝酸 盐 
还 原 为 饺 (DNRA ) 过 程 是 NH; 产生 的 重要 机 制 ,但 
本 研究 中 地 下 水 NH; 浓度 升 高 过 程 伴 随 着 DO 浓度 
升 高 ,并 非 高 度 还 原 环境 ,因此 DNRA 过 程 可 能 不 是 
该 研究 区 NH; 产生 的 主要 因素 ”|。 

8 月 8 日 的 降雨 量 虽 然 较 大 但 并 未 导致 水 力 梯 
度 显著 增 大 ,甚至 大 部 分 监测 井 的 水 力 梯 度 为 负 ， 
呈现 地 下 水 向 河水 补给 , NH; 浓度 也 仅 有 小 幅 上 
升 。 当 降雨 结束 后 ,河水 与 地 下 水 的 水 力 梯度 在 0 
附近 波动 , 即 河水 与 地 下 水 交换 微弱 。 缺 乏 DO 和 
DOC 进 入 河岸 潜流 带 的 情况 不 利于 矿 化 作用 , NH; 
浓度 也 缓慢 下 降 。 表 明 降 雨 事件 通过 改变 河水 -地 
下 水 的 交换 作用 ,控制 营养 物质 的 输入 ,从 而 影响 
河岸 潜流 带 NH; 浓度 变化 。 
3.3 牧区 河岸 带 对 氮 素 的 源 汇 效 应 

夏季 汛期 降水 事件 的 频 发 导致 了 河流 及 地 下 
水 中 氮 素 浓度 大 幅 升 高 ,主要 是 由 于 放牧 导致 大 量 
的 所 负荷 输入 。 研 究 指出 在 放牧 过 程 中 ,动物 摄 入 
70%~90% 的 氮 会 通过 粹 便 和 尿 液 返 回 牧 场 '。 本 
研究 区 牧场 以 牛 群 为 主 , 而 牛 群 放牧 情况 下 的 氮 淋 
失 量 要 高 于 其 他 家 畜 吕 。 此 外 ,研究 区 土壤 以 砂 土 
为 主 , 因 其 孔 际 较 大 渗透 性 较 强 ,更 容易 导致 NO; 
淋 溶 1。 在 整个 监测 期 ,河岸 潜流 带 均 表现 为 NH; 
的 产生 ,而 NO; 的 产生 和 去 除 交替 出 现 (图 8), 因 此 
该 处 主要 考虑 NO; 的 去 除 效 果 。 根 据 公 式 (7) 计 算 
得 到 ,7 月 12 日 (小 十 后 ) 和 8 月 12 日 (暴雨 后 ) 河 岸 
洪流 带 表 现 为 NO; 的 产生 ,是 河流 氮 素 的 源 。 小 雨 
后 左岸 (1L~3L) NO; 的 产生 较 右岸 (1R~3R) 更 明 
显 , 可 能 是 因为 该 时 段 左 岸 DO 相对 右岸 较 高 ,不 
利于 NO; 的 去 除 ( 图 7)。 非 降雨 期 河岸 潜流 带 主 
要 表现 为 NO; 的 去 除 ,去 除 率 为 14%~86% ,是 河 
流 氮 素 的 汇 ,可 能 因为 随 着 落水 期 DO 的 消耗 , 河 
岸 潜流 带 总 体 趋向 厌 氧 环境 (DO<2 mg-L") ,利于 
NO; 的 去 除 '。 并 且 根 据 方差 分 析 结 果 (P=0.25> 
0.05) ,河岸 潜流 带 左右 岸 NO; 去 除 效果 较为 接近 ， 
不 具有 显著 性 差异 (图 8)。 

河岸 带 植 被 群落 物种 较为 丰富 ,可 以 有 效 地 固 
定 和 吸收 多 余 的 氮 素 养分 ,从 而 减少 了 营养 物质 的 
流动 性 ,大 大 增加 了 营养 物质 的 停留 时 间 ,降低 NO; 
WEBS ABS SU. Zuazo 等 3 进行 草地 植被 恢复 后 与 
没有 植被 的 裸露 土壤 相 比 ,发 现 地 下 水 NO; 浓度 减 
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Fig. 8 Production and removal of NO; in riparian hyporheic zones (a) and variance analysis of NO; 


removal level on both banks (b) during the monitoring period 


少 了 35%。 因 此 保护 河岸 带 生 物 多 样 性 ,防止 水 土 
流失 ,是 牧区 毛 素 消减 的 有 效 途 径 。 此 外 ,放牧 不 
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Response mechanisms of nitrate and ammonia nitrogen concentrations to 


hydrological processes in the riparian hyporheic zone of pastoral areas 


XUE Dongyuan, HU Haizhu" ZHANG Jinning, REN Jiawei’ 
(1. School of Ecology and Environment, Inner Mongolia Key Laboratory of River and Lake Ecology, Inner Mongolia 
University, Hohhot 010020, Inner Mongolia, China; 2. Key Laboratory of Ecology and Resource Use of the 
Mongolian Plateau, Ministry of Education, Hohhot 010020, Inner Mongolia, China) 


Abstract: In pastoral and irrigated agricultural areas, nitrogen-containing livestock, poultry manure, and nitrogen 


fertilizers can enter the surface water and groundwater from the soil, and this is the main source of non-point 
source pollution in basins. The riparian hyporheic zone acts as an effective barrier to reduce the nitrogen pollution 
load. Understanding the mechanisms of the migration, transformation, and removal of nitrogen in riparian 
hyporheic zones is key to controlling nitrogen pollution in the whole basin. In this study, an upper reach of the 
Xilin River, located in typical pastoral areas, was selected and its water levels, ammonia ( NH;) and nitrate 
( NO; ) concentrations, as well as the related environmental factors of the river water and riparian groundwater 
during the summer flood season, were continuously monitored. Based on the high- solution measurements, a 
water flow and nitrogen reactive transport model of the riparian hyporheic zones was established using FEFLOW. 
The model fitted using the measured data was found to accurately reproduce the water level dynamics and two 
main nitrogen concentrations in the riparian hyporheic zone. The results indicate that there is a high risk of 
nitrogen pollution in the riparian zones during the summer flood season. The NH; concentration in the riparian 
zones can increase from 0.2 mg + L ' before rainfall events to 7.23 mg + L`' after rainfall events, and the NO; 
concentration can increase from 1 mg: L ' to 8.27 mg-L '. Both measured and simulated results show that the 
nitrogen dynamics in the hyporheic zone are closely related to hydrological processes such as rainfall events and 
groundwater-surface water exchange. During rainfall events, NO; with high mobility was found to infiltrate from 
the river and the ground surface into the riparian zone due to the leaching effect, resulting in a significant increase 
in the concentration. Meanwhile, the groundwater-river water exchange enhanced by rainfall events can further 
regulate NO; and NH; concentrations in the riparian hyporheic zone by controlling the input of nutrients and 
affecting the biogeochemical nitrogen cycles. This study preliminarily reveals the buffering mechanisms of 
pastoral riparian zones in the hydrological and biogeochemical processes involving nitrogen and provides 
scientific references for the nitrogen pollution control in pastoral areas. 

Keywords: riparian hyporheic zone; nitrogen transport; hydrological process; groundwater simulation; 

Xilin River 


